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Abstract 

Employing the Anderson impurity model, we study tunneling properties 
through an ideal quantum dot near the conductance minima. Considering 
the Coulomb blockade and the quantum confinement on an equal footing, 
we have obtained current contributions from various types of tunneling pro- 
cesses; inelastic cotunneling, elastic cotunneling, and resonant tunneling of 
thermally activated electrons. We have found that the inelastic cotunneling 
is suppressed in the quantum confinement limit, and thus the conductance 
near its minima is determined by the elastic cotunneling at low temperature 
{ksT <C r, V: dot-reservoir coupling constant), or by the resonant tunneling 

of single electrons at high temperature (fc^T 3> T). 
PACS numbers: 73.20.Dx, 73.40.Gk 
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During last decade, there has been a rapid advance in the field of single electronics, 
and accordingly much scientific attention has been given to transport properties through 
ultra-small tunnel junctions such as GaAs quantum dot [0-0]. In a quantum dot with small 
capacitance, "Coulomb blockade" of tunneling occurs for small bias voltage V when the 
charging energy (Coulomb energy U) in the dot is sufficiently large as compared to a thermal 
energy ksT '. It occurs because even a single tunneling event increases the electrostatic 
Coulomb energy of the system considerably. However, even in this regime, a finite current 
can flow via virtual intermediate states arising from the quantum fluctuation of macroscopic 
electric charge in the central electrode of the system. This process in a quantum dot, 
so called cotunneling or macroscopic-quantum-tunneling, was first pointed out by Averin 
and coworkers ||, and is considered as setting a limit to the performance accuracy of the 
single electron transistor. They have shown that the transport near conductance minima 
is dominated by the inelastic cotunneling process involving the creation of an electron-hole 
excitation in the central electrode, and predicted an algebraic variation of the leakage current 
with applied voltage (~ V 3 ) and temperature (~ T 2 ). The theory of cotunneling has been 
derived within the lowest order perturbation when the energy discreteness in a quantum dot 
is not important, i.e. the continuous energy spectrum is assumed in the central electrode. 
The inelastic cotunneling has been observed both in metal-insulator-metal tunnel junctions 
|5||| and in a 2D electron system of GaAs/Ga(Al)As heterostructure JF||§. 

When dealing with an ultra-small quantum dot, the effect of level discreteness (energy 
quantization: A) becomes very important. The effect will be more prominent in semicon- 
ductor systems than in metallic systems, due to much lower electron concentration (p) and 
lower effective electron mass (m*) in semiconductor systems (recall that a free electron ap- 
proximation yields A = g ^ F ^ v = vm*(^S P y/ z ^ or systems). There have been quite a 
few experimental evidences exhibiting coexistence of the charge and energy quantization 
in the tunneling properties ||. Furthermore, it is shown that the level spacing A can be 
even comparable to the Coulomb energy U for a Si-based quantum dot transistor |9[]. For 
such systems, the "quantum confinement" will become significant as much as the Coulomb 
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blockade for the tunneling in the single electron transistor. This kind of electronic confine- 
ment could be realized for an ultra-small dot at relatively high temperature. Illustrating 
typical parameters, a Si dot with diameter of 20nm would have A ~ 25 meV, and then the 
confinement of electron could be realized even at T < O(100)K for U ~ 15 meV ||. 

We address in this paper whether the inelastic cotunneling phenomenon is really a lim- 
iting factor in operating single electron quantum dot devices. For this purpose, we have 
examined tunneling properties of an ideal quantum dot coupled to two reservoirs in terms 
of the Anderson impurity model, where the quantum confinement is important as much as 
the Coulomb blockade. Special attention is focused on the temperature dependence of the 
inelastic cotunneling by treating the Coulomb blockade and the quantum confinement on 
an equal footing. We have found that the characteristic of the tunneling in the quantum 
confined system far from conductance maxima is qualitatively different from the case where 
the level discreteness can be neglected. In this case, the "inelastic" cotunneling is substan- 
tially suppressed at low temperature, while the "elastic" cotunneling or the resonant single 
electron tunneling of thermally activated electrons dominates in the system. Therefore it 
is expected that, in small enough quantum dots, there would be no substantial limitation 
on the performance accuracy in practical devices by a macroscopic quantum tunneling of 
charge. 

We start with the simplest Anderson model Hamiltonian to describe an ideal quan- 



tum dot (labeled by D) weakly coupled to two electron reservoirs (labeled by L and R): 

H. = H.L + 7~Ir + Tin + 7~(-t 
n HR) = £ kcl a c ka (1) 

n a n a i 

a a>a' 

7~Lt= iy ka c\ a d a + h.c.) . 

k,a£L,R 

Here 'Hl-i'Hri'Hd and Tix represent Hamiltonians of the left reservoir, the right reservoir, 
an interacting dot, and tunneling between the dot and reservoirs, respectively. The levels in 
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the dot are labeled by an index a, and U denotes the Coulomb interaction between electrons 
in the dot. The states in the reservoirs with energies e k are coupled to the dot by hopping 
matrix element Vk a - The transition rate of an electron between the level e a and the reservoir 
L(R) is given by 

r2 w (w) = 2tt £ \V ka \ 2 8{u-e k ). (2) 

keL(R) 

This Anderson-type Hamiltonian has been recently employed to describe transport in the 



quantum dot structure JXTJ— |T^] . In the Coulomb blockade and the quantum confinement 
limit, we have 

C/,A> k B T,T L(Rh eV (3) 

where eV corresponds to a potential energy difference coming from the bias voltage across 
two reservoirs. 

The cotunneling refers to a simultaneous tunneling of two electrons through intermediate 
states with an extra electron or hole in a quantum dot. This second order cotunneling process 
is called "inelastic" if there remains an electron- hole excitation after the process (Figl.(a)), 
whereas the process is said to be "elastic" if no electron-hole excitation is left (Figl.(b)). 
The inelastic cotunneling current can be simply calculated by using the Fermi golden rule. 

The initial eigenstate \I) of the Hamiltonian 7Y = + 7~Lr + 'Hd can be written 
as \4>l)\^ n )\4>r) i where \4>l(r)) denotes the Fermi sea of left (right) reservoirs, and \ip N ) 
represents an iV-particle eigenstate of Tin- There are two kinds of final states for inelastic 
cotunneling processes. Those are \F\) (an electron (hole) tunnels from L(R) to R(L)), and 
\F 2 ) (an electron (hole) tunnels from R(L) to L(R)), such that |iq) = Cppdpcft a Cka\I) , I-P2) = 
c' ka d a alcpj3\I) (k G L,p G R). The inelastic cotunneling yields the current as 

P n (V) = e( ll (V)- l2 (V)), (4) 

where 71 and 72 are statistical averages of the transition rates from the initial state {|/)} to 
the final states and {|F 2 )|, respectively, and are expressed as 



Ji = (^E\mT 1F 1 ^H T \F i }\*8(E I -E Fi )) i = l,2, (5) 
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h F , Ei — Ho 

with (■ • •)/ denoting the statistical average over the initial states. 
For a constant = ^l(r)( £ )-> 71,72 are given by 

7i00 = ^ E M 1 - "«)> / «fcfe7(e)(l - f{e')) 

x ( — h — ] 5(e' -ep + e a -e- eV) (6a) 

\e a -e + /j, a e'-ep + ^J 

1%{V) = 7iH0- (6b) 

Here (• • •) denotes statistical average, and /(e) = l/{er e + 1) is the Fermi distribution 
function for the reservoirs. // a (/•*&) represents the charging energy of the virtual intermediate 
state with one extra electron (hole). They are defined by 

pi a = NU- n L , (7a) 
H b = Hn-(N- 1)U, (7b) 

with //l(//r) being the chemical potential of the left (right) reservoir which satisfies the 
relation \i L — ji R = eV. Equation @ describes well the inelastic cotunneling when the 
system is not too close to the conductance maxima. 

Note that, if one assumes a continuous spectra for electrons in the dot (A — ► 0), Eq.(^jD 
becomes equivalent to the result by Averin and Nazarov [|J, since then (rip(l — n a )) = fp(l — 
f a ) with / being the Fermi distribution function. However, in the quantum confinement 
limit, the inelastic cotunneling process described by Eq.fl^) gives rise to drastically modified 
characteristics due to following two reasons. First, for the electrons in the dot, the Fermi 
distribution function cannot be used in the quantum confinement limit (ksT,eV A), 
because the number fluctuation is very small [|T^,|T^]. Second, more important difference 



results from the fact that there is no available excitation for inelastic cotunneling process 
with energy e such that < e < A. Therefore it is expected that the inelastic cotunneling 
rate is much suppressed in the quantum confinement limit. 
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To see the functional form of the inelastic cotunneling current in the quantum confine- 
ment limit, it will be sufficient to consider lowest two levels. It is because contributions from 
other levels would be exponentially small as compared to those of two levels. Let's assign 
the two levels as e± = and 62 = A with N — 1 for the ground state of the dot. Then we 
have (nin 2 ) = 0, and Eq.(|6]) reduces to 

r L r R 



H(V) = -^ I dsf(s)(l-f(e')) 



( n i) ( a \ , „ +-7^—1 5(e'-s + A-eV) (8) 



A-e + n a e' + n b 
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+ (n 2 ) + 7 5(e' — e — A — eV) 

Here (m) = 1/(1 + e~ pA ) and (n 2 ) = e _/3A /(l + e _/3A ). The first(second) term in Eq.(g) 
describes the excitation (relaxation) process from initially ground (excited) state in the dot. 
Further, in the quantum confinement limit, (n 2 ) is exponentially small, and so a contribution 
from the second term can be neglected. Then using (rii) ~ 1, Eq.([8]) becomes 

lx{V)^^ fdef{e)f{A-eV-e)[ 1 + 1 ) . (9) 

h J \A-e + fx a e — A + eV + HbJ 

Since f(e)f(A — eV — e) is peaked around e ~ A/2, Eq. (0) is approximated by 



^^(^A7^ + ^) hMn*-eV-e) (10) 



^ ( W + W) 2 (A - eV) e^e^ A . 

/i ^a + A/2 Hb-A/2) 1 ; 



Notable in Eq. ^ [10]) is that the inelastic cotunneling current in the quantum confinement 
limit has a temperature dependence e _/3A suggesting very much suppression at low temper- 
ature. This is in contrast to Independence Averin and Nazarov have obtained for A — > 
limit |4]]. The suppression of the inelastic cotunneling in the quantum confinement limit 
originates from the reduction of phase space for the electron-hole excitations in the dot. 
Hence, in this limit, other kinds of processes such as elastic cotunneling or single electron 
tunneling of thermally activated electrons is expected to determine conductance minima 

m. 



The contributions from the elastic cotunneling and the single electron tunneling of ther- 
mally activated electrons can be estimated by Landauer-type formula through an interacting 
region ]IS| . For the Anderson Hamiltonian, the conductance takes the following form in the 



G = T T ffc(-%)T,P*(e), (U) 



linear response regime, 

a = . 

h J \ de 

where f = f{e) = l/{e^^ + 1) (p = p L = p R ), and p a {e) = -±ImG a (e) is local 

density of states of electrons with label a. Equation of motion - decoupling method can be 
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used to get an approximate solution of the Green's function G a |1^J19||. In the present study, 



more subtle Kondo effect will not be taken into account fl3]14}| . Then Green's functions can 
be expressed as 

G ^)~JT1f72> G '( £ '~ e -(A + 1 t + .T/2 - (12) 
where r = T L + T R . 

The conductance near its minima shows a different feature depending on the relative size 
of T and k B T. For T ^> k B T, ( — §f ) in Eq. (|TTj) can be approximated by 8(e — p) and the 
elastic cotunneling dominates (Figl.(b)): 

G "^^b + (A +[ /-^ j- (13) 

On the other hand, for Y <C k B T, there appears a limit in which the 'resonant' tunneling 
of thermally activated electrons (Figl.(c)) becomes most important. In this case, pi(s) ~ 
8(e), p 2 (s) ~ 8(e — A — U), and thus we have 

G * f hf (0) - /'(A + U)} . (14) 

In any case, the inelastic cotunneling process has much smaller contribution to the conduc- 
tance than other processes. Present results can be summarized as in Fig. 2, which shows a 
schematic diagram of off-resonance tunneling processes. We have examined tunneling prop- 
erties in the quantum confinement regime (T/A <C 1), in which the elastic cotunneling (A) 



or the resonant tunneling (B) is dominating. Note that the inelastic cotunneling (C), which 
has been dealt with in previous studies becomes significant only in the opposite limit 
(T/A > 1). 

In conclusion, we have analyzed off-resonance tunneling properties of an ideal quantum 
dot coupled to two reservoirs, in terms of the Anderson impurity model. We have found that, 
in the quantum confinement limit, the inelastic cotunneling of electrons is much suppressed, 
which is contrary to the case of previously studied A — > limit. The suppression of the 
inelastic cotunneling current originates from the absence of phase space for the electron-hole 
excitations in the dot with energy e such that < e < A. Hence, near the conductance 
minima, the transport is governed by the elastic cotunneling or by the resonant tunneling of 
thermally activated electrons. The present result strongly suggests that the usually thought 
performance limitation of single electron devices arising from the inelastic cotunneling phe- 
nomena would be considerably reduced in an ultra-small quantum dot at sufficiently low 
temperature. 
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FIGURES 

FIG. 1. Various types of electron tunneling processes near conductance minima; (a) inelastic 
cotunneling, (b) elastic cotunneling, and (c) resonant tunneling of thermally activated electron. 

FIG. 2. Schematic diagram of showing various off-resonance tunneling processes for different 
regimes defined by the ratios T/A and T/A. A, B, and C correspond to regions in which tunneling 
processes of the elastic cotunneling, the resonant tunneling of activated electrons, and the inelastic 
cotunneling are dominating, respectively. The Coulomb blockade limit U 3> kT is being considered 
here. 
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